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METHOD AND APPARATUS FOR TIME- 
DOMAIN EQUALIZATION IN DISCRETE 
MULTITONE TRANSCEIVERS 

FIELD OF INVENTION 

The invention pertains to transceivers and modems. More 
particularly, the invention pertains to asymmetric digital 
subscriber line (ADSL) modems used to achieve very high 
speed data communication via telephone networks. 

BACKGROUND OF THE INVENTION 

There is an ever present desire to maximize the speed of 
digital communications via networks, and particularly tele- 
communication networks, such as public telephone systems. 
Accordingly, telecommunication network providers now 
oflfer to their customers many options for coupling to the 
telephone network in addition to the standard analog based 
connection commonly referred to as POTS (Plain Old Tele- 
phone System). Some of the options that are widely avail- 
able are integrated services digital network (ISDN), T-1 
lines, E-1 lines, digital subscriber lines (DSL) and asym- 
metric digital subscriber lines (ADSL). ADSL's can provide 
very high data speeds such as on the order of several 
megabits per second, over a standard twisted wire pair. 
Unlike the traditional data modems used for analog com- 
munication with a telephone central oflBce via a twisted wire 
pair, ADSL requires modems both at the subscriber end and 
at the telephone company Central OfiSce end. Current ADSL 
systems employ discrete multitone (DMT^ technology to 
implement high bandwidth commumications, such as for 
digital TV broadcast, on demand video, high speed video 
based internet access, work at home digital file transfer, 
teleconferencing, home shopping, and information services 
over existing twisted wire pair telephone Hnes. 

The international telecommunications union (TTU) has 
promulgated a standard for ADSL that is commonly termed 
G.lite and which is set forth in ITU-T specification G.9922, 
incorporated herein by reference. Another standard promul- 
gated by ANSI is commonly termed Heavy ADSL and is set 
forth in ANSI specification T1.413, issue 2, also incorpo- 
rated herein by reference. G.lite has 256 samples with 128 
tones (32 tones for upstream communications), each tone 
having a real and imaginary portion. G.lite uses a cyclic 
prefix of length 16 samples. Heavy ADSL has 512 samples 
with 256 tones (32 tones for upstream commimications), 
each tone having a real and imaginary portion. Heavy ADSL 
utilizes a cyclic prefix of length 32 samples. 

FIG, lA is a block diagram of the basic ADSL modem 
functions in accordance with ITU-T G.992.2 with 256 
samples per symbol. The upper half of the diagram repre- 
sents functions in the transmit direction while the lower half 
represents functions in the receive direction. 

In the transmit direction, digital data is generated in block 
102, That data is processed through a scrambler 104 and then 
through a forward error correction (FEQ encoder 106 which 
adds syndrome bytes to the data that will be used for error 
correction by the receiver at the receiving terminal. Next, as 
shown in block 108, the transmit data is encoded using 
quadrature amplitude modulation (QAM). The data is then 
converted from the fi-equency domain to the time domain via 
inverse fast fourier transform (IFFT) 110. 

A 1:4 interpolator 112 interpolates the output of IFFT 
block 110 to produce 256 samples from the 64 samples 
output from block 110. Acyclic prefix is added to each frame 
in block 114. The cyclic prefix comprises data added 
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between the symbols to avoid inter symbol interference 
(ISI). The cyclic prefix is of a standardized length, e.g., 16 
samples for ITU-T G922.2 or 32 samples for ANSI-T1.413, 
issue 2, and is removed at the receiver to recreate the original 
5 transmitted data. The data is then forwarded to a coder/ 
decoder (CODEC) 116. The CODEC encodes the data for 
transmission over the twisted wire pair to the receiving 
device. 

In the receiver portion of the transceiver 100, the received 

10 signal is passed from the twisted wire pair through the 
CODEC 116 where it is decoded. It is then passed to a time 
domain equalizer fTEQ) 118 to shorten the chaimel impulse 
response. Then, in 120, the cyclic prefix is removed. Next, 
an echo canceller (EQ 134 creates an echo cancellation 

15 signal based on the transmit signal which is subtracted by 
subtractor 121 from the receive signal in order to cancel any 
echo of the transmit signal that might return over the twisted 
wire pair and interfere with the receive signal. The echo 
compensated signal is converted back to the frequency 

20 domain by fast fourier transform (FFT) in 122. Then, in 124, 
frequency domain equalization (FEQ) is employed to com- 
pensate for the channel distortion and inter symbol interfer- 
ence (ISI). The receive signal is then processed through a 
quadrature amplitude modulation decoder (QAM) 126 to 

25 decode the tone signal into digital data. That is followed by 
forward error conection (FEC) 128 which uses the syn- 
drome bits that are added by the transmit path FEC encoder 
106 to perform forward error correction. Finally, the data is 
descrambled to extract the true data signal in 130 and then 

30 forwarded to the receiver 132. 

FIG. IB is a block diagram of the basic modem fiinaions 
in accordance with ANSI-T1.413, issue 2. Note that the 
functions are similar to those for ITU-T G.992.2. Notable 
differences include that the receive path FFT 122fl is a 512 
(rather than 256) point FFT and the transmit path interpo- 
lator 112a is a 1:8 interpolator (rather than a 1:4 
interpolator). Both differences are a result of the different 
number of samples between the two protocols. 

FIG. 2 is a block diagram which helps illustrate time 

^ domain equalization in DMT in more detail. In particular, 
block 202 encompasses functions performed at the trans- 
mitter. Block 207 encompasses functions that occur in the 
data channel and block 212 encompasses functions per- 
formed at the receiver. Block 206 encompasses those func- 
tions which result in channel shortening. Thus, it can be seen 
that, in the transmitter 202, a DMT signal, x(n), is generated 
for transmission via the channel. The channel 207 has a 
channel impulse response, h(n), represented by block 208. 
Accordingly, the signal, y(n), received at the receiver 212 is 
the convolution of the original signal x(n) and the channel 
response h(n), i.e.; 

y(n)^nyKn) Eq- (1) 

The time domain equalization circuit 118 in the receiver 
55 212 convolves the received signal, y(n), with a coefficient 
l+a(n). The output y'(n) firom the time domain equalization 
domain circuit 118 therefore is: 

>'(«H</.)M+a(n))=«(n)*(A(«)-(l-Hi(n))) Eq. (2) 

60 The TEQ circuit 118 shortens the channel by the length of 
the cyclic prefix, e.g., 16 samples for G.lite or 32 samples for 
Heavy ADSL. Then the cyclic prefix is removed from the 
signal yXn), as shown in block 120 in FIG. 1. If we call h'(o) 
the channel impulse response corresponding to the shortened 

65 channel, then h'(D)=h(n)*(l+a(n)). 

If the impulse response length of the channel after short- 
ening (hereinafter the shortened channel length, M) is less 
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than or equal to the leagth, L, of the cyclic prefix, i.e., L, 
then each y*(n) will have no dependency on any transmitted 
samples before x(n-L). 

After the cyclic prefix is removed, the actual symbol 
(without inter symbol interference) should remain. 
However, if M^L, then the sampling signal y(n) is depen- 
dent on the transmit sampling signal that appears before 
x(n-L) and, thus, the removal of the cyclic prefix would not 
prevent inter symbol interference and thus degrade the 
performance of the transceiver. Hence, it is desirable to 
design an ef&cient TEQ shortening filter, l+a(n), that makes 
M^L and for which the energy in the channel impulse 
response h'(n) is as small as possible for n^L and as large 
as possible for □<L. 

SUMMARY OF THE IKVENTION 

The invention is a method and apparatus for determining 
the coefficients to be used by a time domain equalization 
circuit in a receiver in a communication system using a 
discrete multitone (DMT) protocol. The method and appa- ^0 
ratus is particularly efficient in that the order of the channel 
response equation after channel shortening is less than or 
equal to the allotted cyclic prefix length utilized in the 
communication protocol. Further, the invention includes an 
efficient method and apparatus for determining symbol 
synchronization for the time domain equalization. 

More particularly, the invention is a method and apparatus 
for generating coefEcients for performing time domain 
equalization on a digital multitone signal involving (1) 
synchronizing the received DMT signal to a corresponding 
signal that was transmitted from said transmitter prior to 
inclusion of any channel response, (2) determining time 
domain equalization coefficients via 



where 



R^kyE{y{n)y{n^k)}-Rjrk) 



and 



FIG. 4 is a flowchart illustrating symbol synchronization 
in accordance with the present invention. 

DETAILED DESCRIPTION OF THE 
5 INVENTION 

The goal of the present invention is to design a TEQ 
shortening filter that produces a new overall channel 
response h*(n) for which M^L. 

10 The impulse response h(n) of the channel has the transfer 
function: 



T represents a transpose function; and 
(3) time domain equalizing the received signal using said 
determined coefficients. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. lA is a block diagram illustrating the functions of an 
ADSL modem transceiver in accordance with ITU-T G 
992.2. 

FIG. IB is a block diagram illustrating the functions of an 
ADSL modem transceiver in accordance with ANSI-Tl. 413, 
issue 2. 

FIG. 2 is a block diagram illustrating in more detail 
channel shortening in discrete multitone communications. 

FIG. 3 is a block diagram illustrating the functions of an 
ADSL modem transceiver for digital multitone communi- 
cations in accordance with the present invention. 



60 



B[Z) 



\Li~lY- 



Eq. (3) 



where L is the length of the cyclic prefix, K is the length of 
the TEQ filter having impulse response (l+a(n)) or A(z). In 
accordance with equation 3, the transfer function of the 
shortened channel is: 

/f(r)^(zy/(r)-Ab+^z-^+&2r-V . . . +V-i)2'^''''' Eq. (4) 

Thus, the length, M, of the shortened impulse response shall 
be set equal to the length of the cyclic prefix, L, in order to 
eliminate ISI between two consecutive symbols by remov- 
ing the cyclic prefix to improve the overall performance of 
the transceiver. 
Let 



Eq. (5) 



and 



35 nnH<n)x(n^l)x(n-2) . . . x(/,-(l-l))y(/,-%(n-2) . . . y{n-{K~ 

1)]^ Eq. (6) 

where T denotes the transpose of a matrix. From equations 
(3) through (6), the estimated signal, y(n), at the receiver is: 



and the error is: 



Eq. (7) 



Eq. (8) 



From the well known LMS (Least Mean Square) error 
algorithm, the optimal solution I^^, is 



50 



where 



55 



R'=E{W(n)W^(n)} 



Eq. (9) 

Eq. (10) 
Eq. (11) 



E{ } in the equations above represents a generalized aver- 
aging fiinction or an ensemble average. must be gener- 
ated before a value for a can be found. 
Let 



RJ^k)^E{x(ri)x{n-^k)}^Rj,-k) 
R,^k)^{x{n)y(n4-k)}^^{^k) 
65 Ry^k)«£My{»^k)}^RJ.-k) 
With equations (9) to (14), we have 



Eq. (12) 
Eq. (13) 
Eq. (14) 
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Equation (15) can be rewritten as 



20 



where 



RM 



R:^iL-\) /f«CL-2) 



/f^CL-2) 

«xx(0) 



RyyiO) 



,(2) 
,(1) 



Ryy{K-2) RyyiK-3) 



RxyiO) 



/f«{-2) 



R,yiL-2) /f^(L-3) 



/f^(A"-2) 
/f^(0) 

R,yi-iK-l)) 

R.,(~iK-2)) 
R^iL~K) 



Eq. (16) 



Eq. (17) 



Eq- (18) 



35 



Eq. (19) 



and 



r^-[R.yiO)R,y(l)R,^{2) . . . R^iL~l)T 



Eq. (20) 



Eq. (21) 



In the equations above, and R^^ are positive definite 
symmetric matrices with the dimensions LxL and pC-l)x 
(K-1), respectively. R^, and rj are Lx(K-l), Lxl and 
(K-l)xl matrices, respectively. We are interested in finding 
the coeflficients of A(z), that is a^,. Hence, from equation 
(16), we have 



RJ>.prR.yP.^r^ Eq. (22) 

RxyKprRy^oprrz- Eq. (23) 

From equation (22), we can solve b^, as 

KfM'^jur\r,-^R.^^i^ Eq. (24) 



60 



By substituting equation (24) into equation (23), we have 

{R^^R^-^R,y'Ry,)a^rr2-(^./R^''r;) 



^op^i^xy^^xx ^Rxy-^)yT\'" ^'^x/^xx'^'" l)- 



Eq. (25) 



Equation (25) can be solved in several ways. In order to 
reduce computation time and memory usage, some condi- 
tions will be considered and some matrices can be pre- 
25 calculated and stored in a table. In equation (25), R^ is an 
autocorrelation of the training signal x(n) that is prendefined 
in ITU-T G.992.2 and is a positive definite symmetric 
matrix. Hence, the matrix (R^^R^'^R^-R^^) also is a 
positive definite matrix. The implementation for equation 
3Q (25) can be summarized as follows: 

(A) Since R^^"^ is a positive definite matrix, it can be 
rewritten as R^'^-V^D^^V, where is a lower triangular 
matrix whose main diagonal elements are all 1 's, and D^^^ is 
a diagonal matrix. These matrices can be pre-calculated and 
stored in a table. Only L(L-l)/2 locations are needed to store 
the V matrix and L locations are needed to store the 
matrix. Hence, equation (25) becomes 



40 



{{VR^yYD^VR^-Ryy)a^r^-{{VR,/D^,Vr^ 



Eq. (26) 



(B) The matrix VR^ which appears 3 times in equation 
45 (26) can be computed once and used in aU three instances. 

The matrix F=((VR^)^D^^ VR^-R^ in equation (26) is a 
(K-l)x(K-l) positive definite matrix. 

(C) By employing the Cholesky decomposition method 
(L. R. Rabincr and R. W/ Schafer, "Digital Processing of 

5Q speech Signals", Prentice-Hall, Inc., Englewood Qiffis, N J., 
1978), equation (26) can be solved for a^^ 

Equation (3) is a general formula that covers cases of 
causality than can occur in the real world. Depending on the 
behavior of a particular design, the transfer function may 

55 alternately be written as: 



B{z) ^ 1 + fcjg-' 4-^22'^ + + ^'iL-nZ'^^^* Eq. (3a) 



and the error e(n) that is to be minimized may be written as: 

e(n)'^nys{n) Eq.(8a) 

where x is the estimated signal at the receiver. The same 
65 concepts described above to arrive at Eq. (25) based on Eq. 
(3) can be used with respect to Eq. (3a) to derive a similar 
solution. 
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FIG. 3 is a block diagram illustrating the operation of a in equation (15) should be computed. TXvo major factors 

modem transceiver for digital multitone communications in in determining the symbol synchronization are the length of 

accordance with the present invention. It is essentially the TEQ filter and the characteristics of 0^{ni). The length 

identical to FIG. 1, except for the new TEQ filter 302 ofthe TEQ filter determines how many peaks and valleys of 

implementing the operation described hereinabove. 5 0xy(ni) will be required to get the desired performance of the 

To effectively implement the new computation of TEQ TEQ coeflBcients. The characteristics of 0^(m) determine 

coefScients in accordance with equation (25) or (26). y(n) the spacing between the peaks and valleys, 

from the CODEC should be well synchronized to the A flowchart to implement the symbol synchronization in 

transmitted signal x(n). If the received signal y(n) is not well accordance with the description above is shown in FIG. 4, 

synchronized to the transmitted signal x(n), then R;^(k) in lO In step 402, the impulse response of the channel is 

equation (15) may not generate the correct value and the determined by utilizing FFT techniques or cross correlation 

TEQ coefficients obtained would not be an optimal solution. techniques. In step 404, the impulse response 4);j^m) for 

Accordingly, it is important to accurately establish symbol 0=m^N passes through an energy detector with a window 

synchronization to achieve the best performance of the of the size of the channel prefix length, termed CPLEN, in 

above algorithm. 15 order to determine the maximal energy and the range of the 

A technique for symbol synchronization will now be channel signal. In step 406, the noise floor of the channel is 

discussed. Let us consider handshaking between two DSL determined from the entire chaimel impulse response. The 

modems. The Central Office (CO) modem sends the average noise floor can be calculated easily by avoiding the 

C-REVERB2 signal (see ITU-T G.992.2 specification for signal energy around the maximal energy area, 

definitions of variables) to the subscriber's modem. The 20 In step 408, the start point of the channel signal is 

subscriber's modem uses the C-REVERB2 signal to perform determined by using the noise floor condition. Then in step 

symbol synchronization and compute the TEQ coefiHcients 410, starting at the start point determined in step 408, the 

simultaneously. In order to find the symbol synchronization, values of the first C peaks and valleys that are totally within 

a cross correlation technique is employed as the channel impulse response are determined, where C is a 

25 pre-determined constant based on experiment. For example, 

N ^ Eq. (27) in the implementation of rrU-T G.992.2 for ADSL modems, 

b^im) = YjXin)y{n-^m)= V jc(rt)^ + m - it) CPLEN is 16 and the length of the TEQ is set to 6. With 

»-o ^-^ *-o these conditions, the number of peaks and valleys, C, can be 

^ set to 6. 

Z Y* L 30 In step 412, equation (25) is used to solve for the TEQ 

h{k)^M,nU,n^m-k) coefiSdents starting at the sample that is CPLEN samples 

before the start point determined in step 408. In step 414, the 

^ coefiScients are returned for use in the TEQ filter. 

= 2j *C*)^j"('" - *) While, as an ahemative, it is easy to compute all possible 

35 sets of TEQ coefiScients and then choose the set that pro- 
vides the best performance, such a technique is time con- 
where N is the number of samples in one DMT symbol, suming. On the other hand, the technique described herein 
0^(m) is the cross correlation of x(n) and y(n+m), and for determining the TEQ coefficients is much less burden- 
0^{m) is the autocorrelation of x(n) and x(n+m). Since FFT some. 

(Fast Fourier Transform) is an indispensable module in a 40 The above-disclosed algorithm may be implemented in 

DMTDSLmodem, we can save memory and computational any of the known ways of implementing TEQ filters, 

overhead by calculating ^jji^m), O^m^N, m the frequency including, but not limited to, a programmed microprocessor, 

domain (i.e., hKl). By converting both sides of equation 27 a micro controller, a slate machine, combinatorial logic or 

by FFT, we derive: any combination ofthe aforementioned. Thus, in accordance 

45 with the present invention, an efficient method and apparatus 

<i>x/*i>)-OTtM'»))w/(to)<p„(a)) Eq. (28) for calculating the coefficients for a time domain equaliza- 

u u/ \ • *u TTirr r<L u i un\T*/\ ^on filter is provided along with an efficient frequency 

where H(a)) is the FFT or the charmel response h(k) Let x(n) , . ^. t_ i ^ - 

u *L nrjT/rinTio * i J n/ \ r-r-r/ / \\ domam apparatus and mcthod f Or symbol Synchronization, 

be the C-REVERB2 signal and P(co)=FFT(x(n)). Then equa- it • fi. j ^ c i ^ j- ^ ^ 

/^o\ I. V \ -1 Havmg thus described a few particular embodiments of 

tion (28) can be wntten as: • i. jve j 

^ ^ 50 the invention, various alterations, modifications, and 

<&,^(o))=i/Ca))4>j„((i))=/^(<D)PCa>)p*(<i>)=y[o)V*(o>) Eq. (29) improvements will readily occur to those skilled in the art. 

Such alterations, modifications and improvements as are 

where P*((o) is the conjugate of P((o) and Y(a))-FFT(y(m)) made obvious by this disclosure are intended to be part of 

in equation (29). From equation (29), the cross correlation of this description though not expressly stated herein, and are 

0^(m) for Oim^N can be easily found by employing an 55 intended to be within the spirit and scope of the invention. 

IFFT (Inverse Fast Fourier Transform) module as Accordingly, the foregoing description is by way of example 

<► i/n)=/Fir7T<»«(<o)) E (30) limiting. The invention is limited only as 

° ^ defined in the following claims and equivalents thereto. 

0^(m) for O^m^N in connection with a subscriber's DSL I claim: 

modem also can be interpreted as the impulse response of 60 1. A method of generating coefficients for performing lime 

the whole channel, including (1) the transmitter filter of the domain equalization on a digital multitone signal compris- 

CO DSL modem, (2) the characteristics of the two wires ing a plurahty of symbols received from a transmitter over 

between the CO modem and the subscriber's modem, and a communication channel, said method comprisiag the steps 

(3) the receiving filler in the subscriber's modem. of: 

There are many factors that determine which sample 65 (1) synchronizing said received digital multitone signal to 

between 0 and N should be selected as the start point. a corresponding signal that was transmitted from said 

Therefore, the cross correlation R^ and the auto-correlation transmitter prior to inclusion of any channel response; 
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(2) determining time domain equalization coeflBcients, 
a^^ for a TEQ filter via 



said received digital multitone signal has a transfer function 
of: 



Biz) 



"of» '»^xt *^xy '*>y> \'2 '^jy ^ 

where 

Rj,k)=E{xin)yin+k)}«R^-k) 
R^k)=£{yin)y(n+k)}^^{-k) 

E{ } in the equations above represents a generalized 

averaging function or an ensemble average 

and 



r = E{yin)Win)). 



where 

W(nHx(n) x(n-l) x(n-2) . . . x(ii-(L-l))y(D-l) y(n- 

2)...y(n-(K-l)r 

L is the length of a cyclic prefix 

K is the length of the TEQ filter 

T represents a transpose function; and 
(3) time domain equalizing said received signal using said 

coefficients determined in step (2). 
2. The method as set forth in claim 1 wherein step (2) 
comprises solving for a^, using one of 



where Vn)^^V=R^-\ 

is a lower triangular matrix whose main diagonal 
elements are all l*s, and 
is a diagonal matrix. 

3. The method as set forth in claim 2 wherein step (2) 
comprises the steps of: 

(2.1) calculating each of and D^^ once; 

(2.2) storing and in a memory; and 

(2.3) retrieving and when needed for determining 

4. Tne method as set forth in claim 2 wherein step (2) 
further comprises: 

(2.4) calculating VRj^^ once and using that calculated 
value for each occurrence of VR^^ in the equation. 

5. The method as set forth in claim 2 wherein step (2) 
comprises: 

(2.5) solving the equation for a^^, by employing a 
Cholesky decomposition technique. 

6. The method as set forth in claim 1 wherein said method 
is used in connection with a digital multitone receive archi- 
tecture that utilizes a cyclic prefix of length L and wherein 
step (2) provides a coefficient wherein M^L, where M is the 
impulse response length of a shortening channel utilized in 
said method. 

7. The method as set forth in claim 6 wherein L=M, 

8. The method as set forth in claim 7 wherein step (2) 
comprises assuming that an impulse response of a channel of 



9. The method as set forth in claim 6 wherein said time 
domain equalization is performed in connection with an 
asynchronous digital subscriber line. 

10. The method as set forth in claim 9 wherein said time 
domain equalization is used in connection wi± a commu- 
nications system in accordance with ITU-T G. 992.2 sped- 

15 fications. 

11. The method of claim 1 wherein step (1) comprises: 

(1.1) cross correlating the received signal with the trans- 
mitted signal to generate an impulse response of the 

2Q chaimel; 

(1.2) determining a maximal energy window of said 
impulse response and a range of said channel response; 

(1.3) determining a noise floor of said chaimel impulse 
response; and 

(1.4) determining a start point of said chaimel impulse 
response using said noise floor. 

12. The method as set forth in claim 11 wherein step (1.2) 
comprises passing said impulse response through an energy 

30 detector with a window of size L samples. 

13. The method as set forth in claim 12 wherein step (2) 
is performed starting with a sample of said received signal 
that is L samples before said start point. 

^5 14. The method as set forth in claim 11 further comprising 
the steps of: 

(1.5) converting said received signal to fipequency domain 
before step (1.1); and 

(1.6) converting said channel impulse response from 
frequency domain to time domain after step (1.1). 

15. A method of recovering a signal that was received 
after being transmitted in digital multitone with a cyclic 
prefix, said method comprising the steps of: 
45 (1) time domain equalizing said signal; 

(2) removing said cyclic prefix fi-om said signal; 

(3) frequency domain equalizing said signal; 
wherein step (1) comprises convolving said signal with 

50 l+a(n) wherein a(n) is generated via 



55 

^op^{Rx/Rxx~^Rxy-Ry^~^ ^-Ra/Rxx'^Ti 

where 

60 R^k)-E[x{n)x{»+f^)}-RJ.-k) 
R,/k)^{x(^)yin^k)}^^-k) 
R^k)^E{yin)yin^k)}^R^-k) 



40 



65 



E{ } in the equations above represents a generalized 
averaging function or an ensemble average 
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and 



where 

W(n)=[x(n) x(n-l) x(n-2) . . . x(n-(L-l))y(n-l) 

y(n-2) . . . y(a-(K-l) f 
L is the length of a cyclic prefix 
K is the length of a TEQ filter 
T represents a transpose function. 
16. The method as set forth in claim 15 wherein step (1) 
comprises solving for using one of 



10 



15 



12 



(4.5) converting said signal to frequency domain before 
step (4.1); and 

(4.6) converting said channel impulse response from 
frequency domain to time domain after step (4.1). 

28. A receiver for a digital multitone modem comprising: 
means for synchroniziDg a received digital multitone 
signal to a corresponding signal that was transmitted 
from a transmitter prior to inclusion of any channel 
response; 

means for determining, for said received signal, time 
domain equalization coefiBcients, a^^ for a TEQ filter 
via 



m^YO,jm^-R„)a^rr:,-{{VR,yYD,^Vrd-. 
Where V^D^V=R^-\ 

is a lower triangular matrix whose main diagonal 
elements are all I's, and 
D^^y is a diagonal matrix. 

17. The method as set forth in claim 16 wherein step (1) 
further comprises: 

solving one of said equations for a^, by employing a 
Cholesky decomposition technique. 

18. The method as set forth in claim 15 wherein said 
cyclic prefix is of a length, L, and wherein step (1) shortens 
a channel of said signal by M, where M^L. 

19. The method as set forth in claim 18 wherein L=M. 

20. The method as set forth in claim 19 wherein step (1) 
further comprises assuming that an impulse response of a 
channel of said signal has a transfer function of: 



21. The method as set forth in claim 18 wherein said 
signal is transmitted via an asynchronous digital subscriber 
line. 

22. The method as set forth in claim 21 wherein said 
method is used in connection with a communications system 
in accordance with ITU-T specification G.992.2. 

23. The method as set forth in claim 15 further comprising 
the steps of: 

(4) determining a starting point for calculating a^^. 

24. TTie method of claim 23 wherein step (4) comprises: 

(4.1) cross correlating the received signal with the trans- 
mitted signal to generate an impulse response of the 
channel; 

(4.2) determining a maximal energy window of said 
impulse response and a range of said channel response; 

(4.3) determining a noise floor of said channel impulse 
response; and 

(4.4) determining a start point of said channel impulse 
response using said noise floor. 

25. The method as set forth in claim 24 wherein step (4.2) 
comprises passing said impulse response through an energy 
detector with a window of size L samples. 

26. The method as set forth in claim 25 wherein step (1) 
is performed starting with a sample of said received signal 
that is L samples before said start point. 

27. The method as set forth in claim 24 further comprising 
the steps of: 



20 



25 



30 



where 

Rj.k>£[x(n)xin+k)}=R^-k) 
R,^k)-E[x(n)y{n+k)}^RJ.k) 
Ry^k)^{y(n)yin^k)}=R^-k) 

E{ } in the equations above represents a generalized 

averaging function or an ensemble average 

and 



35 



r = E{yinmn)] 



■a 



40 



50 



where 

W(n)H:x(n) x(n-l) x(n-2) . . . x(n-(L-l))y(n-l) y(n- 

2) . . . y(n-(K-l)f 
L is the length of a cyclic prefix 
K is the length of the TEQ filter 
T represents a transpose function; and 
means for time domain equalizing said received signal 

using said coefficients. 
29. The receiver as set forth in claim 28 wherein said 
means for determining time domain equalization coefficients 
comprises solving for a^,^, using one of 

(fi^^Di..VR,y-R^)a^r^.Qi^^D^,Vr^ 



55 



where V^D^^V-R^-\ 

V'^ is a lower triangular matrix whose main diagonal 

elements are all I's, and 
D^^^, is a diagonal matrix. 

30. The receiver as set forth in claim 29 wherein said 
means for determining time domain equalization coefficients 
solves the equation for a^, by employing a Cholesky 

65 decomposition technique. 

31. The receiver as set forth in claim 28 wherein said 
receiver is used in connection with a digital multitone 



03/31/2004, EAST Version: 1.4.1 



us 6,678, 

13 

receive architecture that utilizes a cyclic prefix of length L 
and wherein said means for determining time domain equal- 
ization coefiBcients provides a coefficient wherein M^L, 
where M is the impulse response length of a shortening 
channel utilized in said method, ^ 

32. The receiver as set forth in claim 31 wherein L=M. 

33. The receiver as set forth in claim 32 wherein said 
means for determining time domain equalization coefiScients 
assumes that an impulse response of a channel of said 
received digital multitone signal has a transfer function of: 

^ ~ A{z) 1 + air ^ + air 2 + ...-•- a^^_(jj-tA:-i) ' 

34. The receiver as set forth in claim 31 wherein said 
receiver is used in connection with an asynchronous digital 
subscriber line. 

35. The receiver as set forth in claim 34 wherein said 
receiver is used in connection with a communications sys- 
tem in accordance with mJ-T specification G.992.2. 

36. The receiver as set forth claim 28 wherein said means 
for synchronizing comprises: 
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means for cross correlating the received signal with the 
transmitted signal to generate an impulse response of 
the channel; 

means for determining a maximal energy window of said 
impulse response and a range of said channel response; 

means for determining a noise floor of said channel 
impulse response; and 

means for determining a start point of said channel 
impulse response using said noise floor. 

37. The receiver as set forth in claim 36 wherein said 
means for determining a maximal energy window comprises 
means for passing said impulse response through an energy 
detector with a window of size L samples. 

38. The receiver as set forth in claim 36 wherein said 
means for synchroniziDg comprises: 

means for converting said received signal to frequency 
domain before cross correlating the received signal 
with the transmitted signal; and 

means for converting said charmel impulse response from 
frequency domain to time domain after cross correlat- 
ing the received signal with the transmitted signal. 
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